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We investigate the finite-temperature corrections to scaling in the three-state square-lattice Potts 
antiferromagnet, close to the critical point at T = 0. Numerical diagonalization of the transfer matrix 
on semi-infinite strips of width L sites, 4 < L < 14, yields finite-size estimates of the corresponding 
scaled gaps, which are extrapolated to L — > oo. Owing to the characteristics of the quantities under 
study, we argue that the natural variable to consider is x = Le~ 213 . For the extrapolated scaled 
gaps we show that square-root corrections, in the variable x, are present, and provide estimates for 
the numerical values of the amplitudes of the first- and second-order correction terms, for both the 
first and second scaled gaps. We also calculate the third scaled gap of the transfer matrix spectrum 
at T = 0, and find an extrapolated value of the decay-of-correlations exponent, r/3 = 2.00(1). This is 
at odds with earlier predictions, to the effect that the third relevant operator in the problem would 
give ??p atagg = 3, corresponding to the staggered polarization. 
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I. INTRODUCTION 

The three-state Potts antiferromagnet on a square lat- 
tice exhibits a second-order phase transition at T = 0, 
with distinctive properties. Among these is the exponen- 
tial divergence of quantities such as the correlation length 
and staggered susceptibility. 

While earlier studies agreed in pointing to a temper- 
ature dependence of the bulk correlation length in the 
form 



^(f3) ~ (3 p eMvF) , 
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different conjectures were advanced for the values of p, v 
and x, mostly on the basis of numerical work. In par- 
ticular, the value of x was variously estimated as 1.3 
(transfer-matrix results analysed by the Roomany- 
Wyld approximant and Monte Carlo work ||); 3/4 
(conformal invariance arguments coupled with an analy- 
sis of the eigenvalue spectrum of the transfer matrix 
and 1 (further Monte Carlo work H). Later studies 
applying crossover arguments to transfer-matrix data 
taken with an external field H, near the critical point 
T = H = 0, gave x = 1.08 ± 0.13. Additional evidence 
compatible with x = 1, and v — 2, was found via exten- 
sive Monte Carlo simulations [Q. In this latter reference 
it was argued that, although p ~ 1 gave the best fits 
to numerical data, such a logarithmic correction to the 
dominant behaviour was difficult to justify on theoretical 
grounds; also, a value of p — could be made to fit the 
data, albeit with poorer quality than for p = 1. 

A substantial step towards fuller understanding of the 
critical properties of the model was given in Ref. 0. 



Through a mapping to the six-vertex model, where the 
most relevant excitations are vortices, the authors were 
able to find that the bulk correlation length diverges as 
above, with the following exact values for the correspond- 
ing parameters: x = 1, v = 2, p = 0. Further, they es- 
tablished the form of the leading corrections to scaling, 
so that 



f 00 ({3)=Ae 2 ' 3 (l + &/? e -" + ---) 
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The value b = — 6.65 ± 0.11 was calculated, upon consid- 
eration of the stiffness constant of a related model where 
non- vortex defects are the main excitations. Similar re- 
sults were derived for the bulk staggered susceptibility. 
Finally, it was shown that the data of Ref. [?| are com- 
patible with the predictions just quoted. The effects pre- 
viously ascribed to logarithmic corrections could be ex- 
plained once the corrections to scaling, in the form and 
sign predicted, were taken into account. The constant in 
Eq. (|) was fitted to A = 0.121(3) ||, close to the earlier 
estimate A ~ 0.11 - 0.12 for p = in Ref. |. 

Data in Ref. [?] were taken for 2.0 < (3 < 6.0 (corre- 
sponding to 5 < £00 < 20, 000), onlxL lattices with 
32 < L < 1536. Thus, in most cases extrapolation proce- 
dures were used to estimate the L — > 00 limiting values 
of the quantities of interest. 

On account of the exponential divergences, the error 
bars associated to extrapolated quantities turned out to 
increase steeply for lower temperatures. For example (see 
Table 4 of Ref. , the estimate of ^ starts with a rel- 
ative error of 1% at (3 = 2.5, which slowly grows to 3% 
at (3 = 5.2 but then reaches 16% at (3 = 5.9, and 33% at 
[3 = 6.0. Therefore, the picture at the high-/? end of the 
fits to theory in Ref. || is less than entirely clear. 

Our main purpose here is to complement the test of 
Ref. ||, by means of transfer-matrix data generated on 
L x 00 strips of the square lattice. Being essentially 
exact results of numerical diagonalization, our data do 
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not suffer from the fluctuations intrinsic to Monte-Carlo 
studies, allowing one to reach arbitrarily large (3, in prin- 
ciple; instead, owing to limitations in the largest strip 
width accessible (we used 4 < L < 14, L even, with pe- 
riodic boundary conditions across), the most important 
potential source of uncertainties is the L — » oo extrap- 
olation. This drawback is somewhat mitigated by the 
rather smooth behaviour of finite-L data sequences, as 
shown below. 



II. STRIP SCALING AND FINITE-SIZE 
CORRECTIONS 

The choice of quantities to investigate is, in part, dic- 
tated by specific features of the strip geometry; here we 
have chosen to calculate the first and second scaled gaps: 



where e A 



lim 



L 



A - A, 



1,2 (3) 



are the (L-dependent) largest eigenvalues of 
the transfer matrix. At the critical point /3 — oo, con- 
formal invariance J|] states that these quantities give the 
respective decay-of-correlation exponents; in the present 
case, the lowest gap i = 1 is related to the staggered 
magnetization, with associated exponent ?7stagg = 1/3, 
while i = 2 gives the uniform magnetization decay, 
?/ u = 4/3 Jl(], Q. The next relevant operator is re- 
lated to the staggered polarization |ll||, and will be 
briefly discussed in connection with scaling of the third 
gap (i = 3), in Section [V. 

For finite j3 one is off criticality, thus the r\i above are 
not to be interpreted as exponents; nevertheless, they 
are quantities whose difference from the bona fide f3 = 
oo exponents is expected to depend on powers of the 
(suitably defined) distance to the critical point. 

According to finite-size scaling jl2| one must have, with 
= L = oo) given by Eq. (§): 
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Since the finite-si^e corrections here usually are of larger 
magnitude than the fimte-temperature ones, we shall only 
take into account the dominant temperature dependence 
°f Coo (/?), that is, we shall write 



Le 



-2/3 



(5) 



On the other hand, the incorporation of the finite-L ef- 
fects will be done phenomenologically, as explained in the 
following. 

At T = (that is, x = 0), very good convergence 
of finite-width estimates (rji(L)) of ?7 s tagg,u towards the 
exact results (rn) is attained by assuming corrections of 
the form: 



r> m - n- 4- a '° 4- bi0 

T}i(L) - 7]i + — -\-jj 



(6) 



Table I: Finite-L and extrapolated estimates of r/stagg, Tfr- 
The latter are the results of equal-weight fits of data for 



L = 6, 8 • ■ • 14 respectively 
tion (Extr. 1) and to Eq. 



to a single-power (L 
(Extr. 2). 



J ) correc- 



L 

4 
6 

8 

10 
12 
14 

Extr. 

Extr. 



Exact 



^stagg 

0.308785582 
0.321556256 
0.326473031 
0.328860921 
0.330193867 
0.331011103 

0.3331(1) 
0.333303(5) 
1/3 



»7u 

1.47544318 
1.39168002 
1.36528410 
1.35352975 
1.34726477 
1.34352745 
1.3324(3) 
1.333347(2) 
4/3 



These so-called 'analytical' corrections, in powers of L~ 2 , 
are expected to occur for any theory on a strip geometry, 
as they are related to the conformal block of the identity 
operator JliJ. They will be the main corrections, pro- 
vided that no other irrelevant operator with a low power 
arises (as is the case for the three-state Potts ferromag- 
net 1 14, |l5| where an L~ 4 / 5 term is present). In order 
to illustrate how Eq. (H) works, and to give readers the 
opportunity to try their own extrapolation procedures, 
Table | gives our finite-L estimates of ?7 s tagg and ?y u , 
together with their respective extrapolations via equal- 
weight least-squares fits of data (we systematically dis- 
card L = 4 data) . Error bars quoted are the standard de- 
viations of the estimated intercepts at L _1 = 0, as given 
by standard least-squares fitting procedures. Before go- 
ing further, it must be stressed that this structure of cor- 
rections to scaling is, in principle, specific to strip geome- 
tries jl3]]; thus it is not surprising that different results 
(namely, corrections to £/L given by BL^ + CL- 5 / 3 ) 
have been found for this same model, also at T = 0, on 
fully finite L x L lattices Q. 

In order to disentangle the finite-temperature correc- 
tions (to bulk behaviour) which are of interest here, we 
shall assume that, for fixed x = L e~ 2/3 one can still write 



T]i(L, x) = rji{x) 



Oi(x) bi(x) 



L 2 



L 4 



(7) 



where lim^^o a>i(x) — aio and similarly for the other IE- 
dependent quantities. In this way we expect to account 
for the explicit L-dependence of our finite- width results, 
being left only with that given through the argument of 
Eq. (||), which is intrinsic to scaled gaps. 

We illustrate the validity of the smoothness assump- 
tion just made, by displaying in Table [H] our data for the 
largest value of x used (see below), x max = 0.04096. 

Comparison with Table | shows that, although stan- 
dard deviations have increased by roughly two orders 
of magnitude, they still keep within quite reasonable 
bounds, giving credence to the smoothness assumption 
underlying Eq. (0) for all intermediate-x values used 
here. 
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Table II: Finite-!/ and extrapolated estimates of 771, 772, for 
x — 0.04096. The latter are the results of equal-weight fits of 
• 14 respectively to a single-power (L~ 2 1 

1) 



data for L = 6, 
correction (Extr. 



• 14 respectively to a single-power (L ) 
and to Eq. (0) (Extr. 2). 



10 
12 
14 

Extr. 

Extr. 



0.395983934 
0.402292849 
0.404971786 
0.406303355 
0.407011253 
0.407389676 
0.4086(1) 
0.40862(6) 



V2 

1.64309174 
1.58471908 
1.55458116 
1.53956140 
1.53059628 
1.52458596 
1.512(1) 
1.5090(6) 



III. FINITE-TEMPERATURE CORRECTIONS 

In the analysis of the extrapolated (bulk) quantities, 
we shall check for corrections to scaling in the x variable, 
that is, 



r)i(x) = r/i(0) + d <Jx + A (V^) 2 + 



(8) 



Note that a literal translation of Eq. (g) would suggest 
that the corrections in Eq. (||) should depend on \fx In x 
rather than y/x alone; however, consistently with the ar- 
gument used in establishing Eq. (|[), here we shall deal 
only with the dominant terms. 

We have taken x-values decreasing by powers of two, 
from x max = 0.04096 to x min = x max /2 27 = 3.05176 x 
10~ 10 . Using, as a first-order approximation, £oo(/?) = 
Ae 2 @ with A = 0.121(3) 13], the above values of x corre- 
spond to the range (3 = 2.29, ^ ~ 10 (x = x maX) L = 4) 
to /3= 12.27, £oo ~ 5.4 x 10 6 (x = x min , L = 14). The 
lower limit was set by determining when the difference be- 
tween the central estimates r?| xt (a; m i n ) — 7y° xt (0) became 
of the same order as the standard deviation of either ex- 



trapolated quantity (see Table |I]] below, where one sees 
that, although this criterion has been followed strictly for 
772, the three smallest-a; entries for 771 are in fact below 
the threshold; however, by performing analyses with and 
without the corresponding data, we have checked that 
this is of no great import to our conclusions) . 

Table II] also shows that, although our extrapolations 
are very precise, owing to the remarkably smooth varia- 
tion of data against L, they seem to suffer from a slight 
lack of accuracy. Indeed, for x — our central estimates 
77| xt stand respectively 6 and 8 standard deviations away 
from the known exact values for 771 and 772. We ascribe 
this effect to systematic errors coming from: (i) the short- 
ness (in L) of our data series, and (ii) higher-order correc- 
tions, ignored in Eq. ((?]). Since, at least for x = 0, such 
errors amount to small differences in the central estimates 
(respectively —0.01% and +0.001% for 771 and 772) rela- 
tive to exact values, and assuming this scenario to carry 
over, continuously and smoothly, to x ^ 0, we shall do 
as follows. In Eq. (|J), for instance, we shall use ??^ xt (0) 
instead of the exact ??i(0); this way we expect systematic 
errors to cancel to a large extent, when considering the 



Table III: Extrapolated values of 771(2;), 7/2(2;). For each x 
they are the result of an equal-weight fit of data for L — 6, 
8---14toEq. (0). 



x 
0. 

3.05176E-10 
6.10352E-10 
1.2207E-09 
2.44141E-09 
4.88281E-09 
9.76563E-09 
1.95313E-08 
3.90625E-08 
7.8125E-08 
1.5625E-07 
3.125E-07 
6.25E-07 
1.25E-06 
2.5E-06 
5.E-06 

1. E-05 

2. E-05 
4.E-05 
8.E-05 
.00016 
.00032 
.00064 
.00128 
.00256 
.00512 
.01024 
.02048 
.04096 



.333303(5) 
.333305(5) 
.333306(5) 
.333306(5) 
.333308(5) 
.333309(5) 
.333311(5) 
.333315(5) 
.333319(5) 
.333326(5) 
.333335(5) 
.333347(5) 
.333366(4) 
.333392(4) 
.333429(4) 
.333483(2) 
.333560(1) 
.333673(2) 
.333839(4) 
.334087(8) 
.334465(13) 
.335054(21) 
.335999(32) 
.33756(5) 
.34025(6) 
.34507(9) 
.35409(11) 
.37183(12) 
.40862(6) 



„cxt 
»72 

1.333347(2) 
1.333352(3) 
1.333354(3) 
1.333358(3) 
1.333362(3) 
1.333369(3) 
1.333378(3) 
1.333391(3) 
1.333409(4) 
1.333435(4) 
1.333471(5) 
1.333523(6) 
1.333596(7) 
1.333701(10) 
1.333849(13) 
1.334060(17) 
1.334363(22) 
1.334796(31) 
1.33542(4) 
1.33634(6) 
1.33768(8) 
1.33969(11) 
1.34276(15) 
1.34754(22) 
1.35523(30) 
1.3681(4) 
1.3904(5) 
1.4309(6) 
1.5090(6) 



difference rj^ xt (x) — 77™* (0). 

Our first test is a single-power fit to scaling corrections: 
we assume 



mix)- 77,(0) =CiX u 



(9) 



and vary u within a reasonably broad range, checking the 
behaviour of the x 2 of the corresponding least-squares fit. 
Our results, using as input the upper half of Table QI 
(14 data plus the x = line, for each fit) are displayed in 
Figure |l|, where very sharp minima can be seen slightly 
above x — 0.5 (to three decimal places, they are located 
respectively at x = 0.508 for 771, x = 0.503 for 772)- This 
signals that (i) corrections depending on \fx are definitely 
present, thus supporting the assumption made in Eq. (||); 
and (ii) higher-order terms are not negligible. Indeed, 
inclusion of data for larger x causes the \ 2 to increase 
steeply, while the sharpness of the dips deteriorates, and 
their position shifts towards larger u. 

Having ensured that square-root corrections to scaling 
are an essential element of the picture, we attempt to 
include higher-order terms, in the manner of Eq. (||). We 
plot (rji(x) — r/i(0))/y/x against y/x, thus one expects: 



= Ci+DiVE + 0(x) 



(10) 
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Figure 1: Semilogarithmic plot of x against u for least- 
squares fits of Eq.(|^). Crosses: r/i; squares: r/2. Data for 
r/2 shifted upwards by a factor of 100 on vertical scale, to 
avoid superposition. For each fit, only data in the interval 
3.05176S - 10 < x < 2.5E - 06 of Table |n] were used (see 
text). 



Table IV: Estimates of amplitudes d, Di (see Eqs. (|8j) -(|1__ 
and x 2 P er degree of freedom (x 2 /DOF) for respective fits 

d Di x 2 /DOF 

Eq. (§), i=l, u=0.5 0.079(4) - 1.5 x 10"' 

En. ©, i=l 0.077(1) 1.18(3) 2 x 10~ 2 

Eq. (§), i=2, u=0.5 0.317(5) - 3.0 x 10~ ; 

Eq. @, i=2 0.313(3) 2.3(1) 2 x 10" 3 



Figure 2: Plots of (r]i(x) — r]i(0))/i/x against *Jx, for i — 1 
(a) and i = 2 (b). Straight lines are linear least-squares fits 
to the subset of data in plot, corresponding to 5.E — 06 < x < 
0.00256 in Table 



III (see text). 



quote only results for u = 0.5. Although, as explained 
above, these do not correspond to the respective absolute 
minima of x 2 , they exhibit a very good quality of fit, and 
it seems more appropriate to compare them (instead of 
those obtained at minimal x 2 ) to the estimates of Cj from 
Eq. (|i"(i|), where the power 1/2 is fixed from the start. 



IV. SCALING OF THIRD GAP 



and attempts straight-line fits. Results are in Figure [|. 
The subset of data considered now is complementary to 
that used in the earlier single-power fits, as higher-order 
terms become more important for x not very small. We 
noticed that inclusion of data from the last 3-4 lines of 
Table |l| caused a quick deterioration of the quality of 
linear fits (the resulting curvature can be seen by naked 
eye); this is probalbly the effect of third- and higher- 
order terms in y/x. An alternative source of errors would 
be the multiplicative logarithmic terms, mentioned in 
connection with Eq. (||) above, and not considered in 
the present approach. Therefore, we decided to keep 
to the range of x for which good linear plots were ob- 
tainable, while using as many data as possible (in order 
to reduce the spread in the estimates of Cj and Di of 
Eq- (Jf))- The best compromise was found by taking data 
for 5.E — 06 < x < 0.00256, shown in Figure || together 
with the corresponding least-squares fits. 

Our estimates of the amplitudes Cj and Di are shown 
in Table IV, where for single-power fits of Eq. (EJ) we 



Finally, we have investigated the scaling of the third 
gap at T — 0. According to theory Jl(], fj"l|| , at the critical 
point there are only three relevant operators, correspond- 
ing (in decreasing order of relevance) to staggered magne- 
tization, uniform magnetization and staggered polariza- 
tion. Although, as recalled above, there is widespread 
agreement between theory and numerical work as re- 
gards the first two, the prediction of Ref. |ll], namely 
that the corresponding decay-of-correlations exponent is 
^Patagg = 3, appears not to have been numerically tested 
so far. ( In Ref. [ll], Monte Carlo simulations were per- 
formed for the respective susceptibility, which according 
to the scaling law ^/v — 2 — r\ is expected to approach 
a constant, with corrections oc L~ A , A = 1, if T] = 3; 
the approach to a constant was indeed verified, while the 
best fit was for A ~ 0.75 instead of unity). 

We have calculated descending eigenvalues of the 
transfer matrix; it would seem plausible to associate the 
third scaled gap to the staggered polarization, especially 
since only three relevant operators are expected to come 
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Table V: Finite-!/ and extrapolated estimates of 773- The 
latter is the result of an equal-weight fit of data for L — 6, 
8--- 14 respectively to a single-power (L~ u ) correction; the 
corresponding u — 1.8 was chosen to minimize the x 2 (see 
text). 

L 773 
6 1.74149553 
8 1.84852826 
10 1.90108612 
12 1.93052844 
14 1.94860339 
Extr. 2.00(1) 

up, and the relationship of the other two to the first two 
gaps is well-established. In order to check self-consistency 
of our results, we used both a standard power-method al- 
gorithm, coupled with Gram-Schmidt orthogonalization, 
and a Lanczos scheme. While for small L and shallow lev- 
els (corresponding to eigenvalues Aj, i < 3) both methods 
gave the same estimates, the Lanczos results displayed 
instabilities for deeper levels and L > 8. At present we 
are not able to explain such discrepancies. Therefore, we 
restrict ourselves to the analysis of the third gap. 

Our results, again displayed in the form 773 = L (Aq — 
A 3 )/7r, are shown in Table [v|. In order to gain an un- 
biased perspective both of the limiting bulk value of 773 
and of the scaling corrections, we attempted a single- 
power extrapolation, 773 (L) = 7/3 + a^/L u with a vari- 
able power u, and monitored the variations of the % 2 of 
the corresponding fits against u. The result was qualita- 
tively very similar to that displayed for the fits of Eq. (||) 
in Fig. [j} a rather sharp minimum, located at u = 1.8 in 
this case, which gave an extrapolated 773 = 2.00(1) (see 
Table [v[ the error bar was calculated by considering the 
estimates on either side of u min = 1.8, for which the x 2 
becomes one order of magnitude larger than at the mini- 
mum). Fixing u=l, inspired by the prediction of Ref. |ll| 
for the susceptibility, gave 773 ~ 2.11. Two-power fits a la 
Eq. (0), using either L~ 2 and L~ A or L^ 1 and L~ 2 also 
gave values between 1.99 and 2.01. There seems to be no 
straightforward way to extrapolate the data of Table |y| 
to include 773 = 3 . At this point we do not know how to 



reconcile our results to the predictions of Ref. |Tl|. 



V. CONCLUSIONS 

In summary, we have undertaken a finite-size approach 
to investigate the corrections to scaling in the three- 
state square-lattice Potts antiferromagnet. Owing to the 
characteristics of the quantities under study, we argued 
that the natural variable to consider is x = Le~ 2/3 . We 
showed that the less-relevant finite-size corrections could 
be accounted for in a phenomenological scheme, based on 
the zero-temperature picture; for the extrapolated scaled 
gaps we supplied convincing evidence that square-root 
corrections, in the variable x, are present, and provided 
estimates for the numerical values of the amplitudes of 
the first- and second-order correction terms, for both the 
first and second scaled gaps. It would be interesting if 
predictions based on theory could be derived, to be com- 
pared with the numerical values of amplitudes obtained 
in this work. 

We have also investigated the behaviour of the third 
scaled gap of the transfer matrix spectrum, and found 
an extrapolated value for the decay-of-correlations expo- 
nent 773 = 2.00(1). This seems incompatible with earlier 
predictions, to the effect that the third relevant operator 
in the problem would give 77p Btagg = 3, corresponding to 
the staggered polarization. 
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